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We study the dimensional crossover from ID to 2D in half-filled and lightly doped, weakly 
interacting A^-leg Hubbard ladders. In this case, the Hubbard ladders are equivalent to a A^'-band 
model. Using renormalization group techniques, we find, that the half-filled ladders exhibit (in the 
spin-sector) an odd-even effect only below a crossover energy Ec oc exp[— a exp(7A*')] (a ^ 1 and 
7 ~ 1 depend on the interaction strength and on the hopping matrix elements): Below Ec, the 
dominant interactions take place within band pairs (j, N — j) [and within the band (A'^ -I- l)/2 
for A'' odd], such that even-leg ladders are an insulating spin-liquid, while odd-leg ladders have one 
gapless spinon-mode. In contrast, above the energy-scale Ec, all bands are interacting with each 
other and the system is a 2D-like (insulating) antiferromagnet; we obtain an analytical expression 
for the Hamiltonian which is similar to the 2D Heisenberg antiferromagnet. Bosonization techniques 
show, that in the charge-sector the Mott insulator is as well below and above Ec of the same type 
as in A''/2 half-filled two-leg ladders. Doping away from half-filling, we find that the effect of an 
increasing doping 5 is very similar to decreasing the number of legs A^: In both cases interactions 
between unpaired bands are suppressed and thereby the antiferromagnetic correlations reduced. 
The resulting band pairs form then insulating spin-liquids and when doped, there is a spin-gap, 
but phase coherence exists only within the band pairs. At higher doping levels 5c ~ 5c{N), phase 
coherence between all band pairs sets in and the system becomes a 2D-like d-wave superconductor 
[Sc{N) ^ for Af ^ oo]. 

PACS numbers; 71.10.Pm, 74.20.Mn 



I. INTRODUCTION 

The physical properties of (quasi-)one-dimensional 
electron systems are often very different from their 2D 
counterpart, for example, in (spin-1/2) Heisenberg an- 
tiferromagnets (AFM): The (quasi-lD) ladders exhibit 
in the groundstate an odd-even effect, i.e., odd-leg lad- 
ders have one gapless spinon-mode, while even-leg lad- 
ders have a spin-gap . The 2D Heisenberg model how- 
ever, has two gapless magnon-modes and long-range or- 
der H,^. Similarly, (quasi-) ID conductors are generally 
Luttinger liquids (LL) while 2D and 3D electron 

systems are (usually) Fermi liquids (FL). 

For the AFM Heisenberg model, the dimensional 
crossover from ID to 2D has been studied in Ref. Q. Via 
the correspondence between number of legs N and tem- 
perature T, 1/T cx N, the nonlinear sigma-model allows 
one in this case to calculate various physical properties; in 
particular, the spin-gap present in even-leg ladders van- 
ishes exponentially as a function of the number N of cou- 
pled spin-chains. Coupled LLs and the dimensional 
crossover to 2D have been investigated by many different 
authors 1-0. One of the reasons was the search for a 
non-FL (i.e., possible LL-like) phase in 2D How- 
ever, no clear evidence for such a phase has been found 
so far 

In the half-filled and lightly doped Hubbard model, 
there is a strong competition between Cooper, umklapp, 
and AFM processes. For an unbiased determination of 
the dominant phase, it is therefore advantageous to use 
the renormalization group (RG) technique In the 



(weakly interacting) 2D case, such RG studies found 
AFM at and close to half-filling and d-wave superconduc- 
tivity for sufficient hole doping [p^[-p7|. The dimensional 
crossover from ID to 2D was only investigated away from 
half- filling (without umklapp and AFM processes) [ p^ . 
The authors concluded that the system becomes a FL for 
TV ^ oo. 

In a previous work, we have studied in detail the 
groundstate properties of the half- filled and lightly doped 
A^ = 3 and 4-leg Hubbard ladders in the case of weak in- 
teractions In contrast to the large U case, where the 
half-filled Hubbard model converges onto the Heisenberg 
model of spin-1/2 |20|, in the small U case, the charge 
degrees of freedom are still present, which gives insight 
in the effect of doping away from half-filling. 

Here, we analyze the crossover from ID to 2D for the 
half-filled (and lightly doped) ladders. We will show that 
depending on the energy-scale (respectively length-scale) , 
the physical properties are very different and that the 
phases associated with these properties are separated by 
crossover energies which are a function of N. In par- 
ticular, the quasi-lD analog of the 2D phase is present 
at higher energies, implying that the crossover between 
different phases can be studied for N finite. An advan- 
tage in quasi- ID compared to the 2D case is the limited 
number of different interactions. This allows us to obtain 
(some) analytical results and to disentangle the interplay 
between AFMs (two gapless magnon-modes), ISLs (insu- 
lating spin- liquid, i.e., only short-range correlations), and 
superconductivity. The 2D limit [N = oo) can then be 
taken within the same phase. 
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For weak interactions, the Hubbard ladders are equiv- 
alent to a A^-band model, characterized by Fermi veloci- 
ties Vj . At half-filling, due to nesting of the Fermi surface 
(FS) , the Fermi velocities of the bands j and N + 1 — j 
are equal (to be called "band pairs" below), such that 
AFM and umklapp processes take place. Using RG tech- 
niques, we find that the half-filled ladders exhibit (in the 
spin-sector) an odd-even effect only below a crossover 
energy Ec oc exp[— a exp(7A'')] (a ^ 1 and 7 ~ 1), but 
are 2D-like (insulating) AFMs above the scale Ec- The 
reason is, that below Ec only interactions within band 
pairs (j, N+l-j) [and within the band (iV -I- l)/2 for N 
odd] are present, while above Ec all bands are interacting 
with each other. We obtain an analytical expression for 
the Hamiltonian which is similar to the 2D Heisenberg 
AFM. Analyzing the charge-sector by bosonization tech- 
niques [^212^, we find, that the Mott insulator is of 
the same type as in N/2 half-filled two- leg ladders. The 
scale Ec decreases fast to zero as the number of chains N 
goes to infinity. The 2D limit is thus an insulating AFM, 
in agreement with what has been found in studies of the 
weakly interacting 2D Hubbard model lI^-|lT 
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FIG. 1. Schematic phase diagram for the (large) A'^-leg 
Hubbard ladders. At half-filling, the system below the en- 
ergy-scale Ec = Ec{N) is an insulating spin-liquid (ISL) and 
above Ec an insulating, 2D-like, commensurate antiferromag- 
net (CAFM). Upon doping away from half- filling, the lad- 
ders become below the chemical potential fi (measured from 
its value at half-filling) a conducting spin-liquid (CSL), i.e., 
phase coherence is only present within pairs of Fermi mo- 
menta {kF,TT — kp) and (tt — kFjkp); above fi, the antifer- 
romagnetic correlations become incommensurate (ICAFM). 
Upon increasing doping, phase coherence between all Fermi 
momenta sets in and the system becomes a 2D-like d-wave 
superconductor. 

Doping away from half-filling (at energies above Ec), 
we find that the effect of a doping ^ > is the same as of 
a finite N < 00: Doped holes suppress the interactions 
between unpaired bands and thereby reduce the antifer- 
romagnetic correlations. The effect of the doping on the 
remaining interactions within a band pair (j, N + I — j) 
is then the same as when doping A'^/2 uncoupled ISLs; 



there is a spin-gap and d-wave-like phase coherence be- 
tween bands j and N + 1 — j. At a higher doping value 
Sc — 5c{N), phase coherence between all band pairs sets 
in and the iV-leg ladders become a 2D-like d-wave super- 
conductor [Sc{N) for TV ^ cxd]. The low-energy 
physics is given by a c?-wave ECS Hamiltonian. The 
origin of the superconducting instability is the same as 
for the 2D case p^-|l^: A Kohn-Luttinger-type attrac- 
tion is generated by AFM processes Further- 
more, while at half-filling the antiferromagnetic correla- 
tions are commensurate, upon doping, they become in- 
commensurate. Fig. |l| gives a schematic overview. 

Summarizing, the iV-leg Hubbard model exhibits be- 
low the length-scale « l/£'c almost the same physical 
properties as the 2D model (i.e., the same correlation 
functions dominate), but in contrast to the 2D case, this 
quasi- ID analog can be treated (partially) analytically. 

We do all our calculations for even N and only refer to 
the odd N case where it is of interest. In Sec. H, we intro- 
duce the TV-leg ladder Hamiltonian, discuss the various 
interactions, and specify along which "path" the dimen- 
sional crossover is studied. In Sec. Ill, we briefly revisit 
the results obtained previously in Ref. p9[ | and analyze 
the half-filled Hamiltonian as a function of the energy- 
scale. We show that there exist for finite N two clearly 
different phases, separated by a crossover energy Ec- We 
derive an (analytical) expression for Ec- In Sec. IV, the 
phase above the energy Ec is investigated in detail and 
shown to be a 2D-like AFM. Finally, in Sec. V, we treat 
the effect of doping away from half-filling. 



II. N-LEG HUBBARD MODEL 

The iV-leg (TV-chain) Hubbard model is given hy H = 
Hq + Tfint, where the kinetic energy is 

Ho = ~tJ2 4six + ^)d^six) + H.c. 

x.,i,s 

-U^4+is(a;)rf«s(a^)+H.c. (1) 
and the interaction term is 



Hint = U 4t (a^)4i {x)dii (x) . 



(2) 



The hopping matrix elements along- and perpendicular 
to the chains are denoted by t and t± respectively and 
[/ > is the on-site repulsion. For our approach, it is 
advantageous first to diagonalize Hq, 



(3) 



where and Vl/js are the creation- and annihilation 
operators for the band j. For open boundary conditions 
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perpendicular to the legs, the dispersion relations ej are 
given by (j = 1, . . . , N) 



ej{k) = -2t cos(fc) - 2t±_ cos [7rj/(7V + 1)] . 



(4) 



The Fermi momenta in each band kpj are determined 
by the chemical potential = £j{kpj) and the filling 
n = 2(7riV)^^ ^Fj- Since we are only interested in the 
low-energy physics, we linearize the dispersion relation 
at the FS, resulting in Fermi velocities Vj = 2tsm{kpj). 
We introduce operators "i! r/ljs for right and left movers 
at the FS. At half-filling, n = 0, and 



2Jt^-tl cos[7:j/{N + l)Y 



(5) 



where j — N + 1 — j . For t = t±, vi = vn ^ 2TTt/N, 
which leads to a singular behavior for large N (the quasi- 
ID analog of the van Hove singularities in 2D). For the 
following, we take tj^ < t in order to avoid singularities. 
Note that the FS is then not flat |2|], but that we still 
have nesting, i.e., kpj + kpj = tt. 

The crucial difference between quasi-lD and 2D are the 
interactions which control the low-energy physics. The 
system is quasi-lD and only a finite number of different 
interactions play a role, provided the energy difference 
between two neighboring bands is larger than the largest 
energy-scale of the system, i.e., te~*'/^ . Therefore, we in- 
vestigate the crossover between different phases (ISL and 
AFM) in quasi-lD, U < IniV, In A^, as a function 
of the energy (this can partially be done analytically). 
The 2D limit. In A^ t/U , is then taken within the same 
phase. In other words, along this particular "path" , the 
dimensional crossover from ID to 2D can be treated in a 
controlled way. 

Since we are interested in the small-C/ low-energy 
physics, we have to take into account only the processes, 
which are present down to arbitrarily small energy-scales. 
For the half-filled case, this implies that only umk- 
lapp processes where the momenta add up exactly to tt, 
kpi -\- kpm = TT are relevant. This condition is only ful- 
filled for the band pairs (j, j), kpj + kpj — tt. When dop- 
ing away from half-filling, the chemical potential (mea- 
sured from its value at half-filling) introduces a low- 
energy cutoff for the various umklapp processes. The 
interactions which have to be taken into account initially 
are then different. 

At half-filling, the interacting part of the Hamiltonian 
consists of forward (f). Cooper (c), and umklapp (u) pro- 
cesses between 2 different bands and (non-)umklapp pro- 
cesses between 4 different bands (in our notation, we fol- 
low Rcfs. [plpl), iJint = H"^^ + H^^, where 



^2B J dx {f[jJmiJLjj - fij^RM ■ Jljj 



+CijJmjJLij 
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+<.jA/^i^)+H-c.] (6) 
where the U(l) and SU(2) currents are defined as 

Jh^J^Y.^i^s'^hJS (7) 

s 

and 

•^'.. = |E*L^-'*/>.v, (8) 



where again h ~ R/L and the are the Pauli matri- 

I ji 

^.(7 n.rr rm i i 

C7. 



ces (p — x,y,z). Due to symmetry, /j^'"^ 
— Cjf. The umklapp currents are 



and 



fil and 
(9) 

(10) 



where e = — ir^. The 4-band interactions involve pro- 
cesses between two different band pairs and {k,k), 



N/2-1 N-j 

j=i k=j+i-' 



(11) 



where 

-^J^ = ^jfefej [jRjkJhk] + JRjkJLkj + H.c.) 
'^"jkkj i-^Rjk ■ ^Lkj + Jfljfe • ^Lkj + H.c.) 



-"jkkj 



jkkj y4,jkhkj + llrkj^Ljk + H.c.) 

jjfefe {^RjJ^Lkk + ^Rkk^LjJ 



H.C. 



(4,.IiS, + i;^_.U,fe+H.c.). (12) 



In our case vj — Vj implying c ' 



p,(7 p,(T 



jkkj jkkj^ jkkj 



jkkf 



Note that the 4-band couplings describe AFM processes, 
i.e., the terms '^^js^Lj 
of TT (along the chains). 



i.e., the terms L]s' result in a momentum transfer 



III. DIMENSIONAL CROSSOVER 

Using the RG method, we study the half-filled A^-leg 
ladders as a function of the energy-scale. We show that 
there exist two clearly different phases, separated by a 
crossover energy Ec <C t. The phase below Ec is an 
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ISL, i.e., 4-band-AFM processes are suppressed and the 
physics is dominated by 2-band umklapp and Cooper 
processes |l9). In contrast, in the phase above Ec (but 
still below the bandwidth t) the 4-band-AFM processes 
are large and dominating (this phase is investigated in de- 
tail in Sec. IV). First (Sec. Ill A), we revisit the results 
for A^-leg ladders obtained previously in Ref. and 
give a criterion for the couplings allowing for a (rather) 
precise distinction between the two phases. In Sec. Ill B, 
we calculate the crossover energy Ec as a function of N. 



(15) 



where v = p,a and /i = ± and the Luttingcr liquid pa- 
rameters (LLP) are 



n 



(16) 



A. ISL phase in even-leg ladders 



In Ref. ||T^, we have shown, that at half-filling (for N 
even) the band pair (j,j) scales towards an ISL and is 
frozen out at the energy-scale te'"^""^^^ , where aj ~ 1. 
The low-energy Hamiltonian is then the same as of N/2 
half-filled two-leg ladders. For increasing N, the initial- 
value U/t (and therefore the energy-scale) , which leads 
to this decoupling (in particular to a suppression of the 4- 
band-AFM processes) and an ISL phase, decreases fast. 
Next, using bosonization techniques ||^,|2l|,^, we show 
that the RG flow of the Luttinger liquid parameter for 
spin triplet pairing, K„j^, allows one to determine the 
energy-scale, where the spin-gap opens. This energy- 
scale can then be interpreted as the crossover energy be- 
tween the ISL phase at low energies and the AFM phase 
at higher energies (which has no spin-gap). 

Bosonization is the method of rewriting fermionic cre- 
ation (annihilation) operators in terms of field opera- 
tors $ct and IIq, satisfying the commutation relation 

a{x) ,\la[y)\ = i5{x — y) (a labels the charge and 
spin-modes for the different bands). The dual field of 
<i>a is defined as dxOa — Hq. The fields contain 
the (physical) charge respectively spin of the particle, 
i.e., the charge/spin-density is pa = dx^a, while the 
fields Oa can be interpreted as the phase of the charge- 
part respectively the spin-part (for ladders, see, e.g., 
Refs. [@|l9[|2|l). We introduce the fields 



(13) 



and 



and 



K, 



<jj± 



/ 47r., + (c- ±/^) 



(17) 



The velocities are u^j^ ~ Vj. The interacting part of 
the TV-leg ladder Hamiltonian contains (apart form other 
terms) the expression [E6| 



(%/4^$„jj+) 



X other cosine terms. 



(18) 



Fields which appear in a cosine become "pinned" in order 



to minimize the energy, e.g., the term — cos(\/47r<I>o-jj+) 
implies (quasiclassicaly) 'I*(Tjj+ ~ and the correspond- 
ing spin-mode acquires a gap. 

The LLPs determine the exponents of correlation func- 
tions |2l|j2^ ] and — together with cosine-terms — the 
pinned fields. Here, as we show in the following, the 
flow of the LLP Kaj+ can be used to determine the 
energy-scale where the spin-gap in the a+ sector opens. 
Note that initially (i.e., at high energies) Kr,j+ > 1 and 
K„j_ = 1 (and Kpj+ < 1, Kpj^ = 1). For the TV-leg lad- 
ders, the RG flow of the couplings is such that w ~t 
and fjj « at low energies, implying Ka-j± < 1. Hence, 
while Ka-j- is always < 1, there is an energy, where Ka-j+ 



becomes smaller than 1, i.e. 
canonical transformation 



■'33 — 



^0-33 



3+VK, 



aj+ 



The 



(19) 



eliminates K^j^ from Hq (it just renormalizes the veloc- 
ities) and results in 



— cos 



^ AirKaj + $ cr j j+ ^ X other cosine terms. (20) 



H 



•^33± 



1 
71 



(H,, ± U^j) 



(14) 



combining the fields of the bands j and ] {v = p, cr), 
e.g., the spin-fields ^aj]± describe triplet/singlet pairs 
between the momenta kpj and kpj- 

The bosonized Hamiltonian takes the following form. 
The noninteracting Hamiltonian Hq plus the contribu- 



tions of the charge-currents /j^ and 
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and the 



component of the spin-currents f^^ and c*^- can be written 



as 



A large Kaj+ therefore tends to depin the field 
while a small Ka-j+ results in a strong pinning of the 
field $o-jj+ (for the dual field ^<Tjj+, it is vice versa). 
The crossover between the two regimes takes place at 
Ka-j+ = 1 (since we are at finite energy, it is a crossover 
and not a transition). The ISL phase therefore requires 
rj < {K„j+ < 1) and rj can be used in the RG flow 
to determine a crossover energy. We note that for the 
pure sine-Gordon model (at zero temperature), the value 
Ka — 1 rigorously separates the massive from the mass- 
less phase [pTt- 
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FIG. 2. The upper figure shows the strong-coupling value 
of the ratio of a 4-band-AFM and a 2-band Cooper cou- 
pling, CI2N -IN / ^In y ^■s a function of the initial-value U/t. 
For increasing A'^, the 4-band-AFM couplings are only sup- 
pressed for very small initial-values. Similarly, the lower fig- 
ure displays the strong-coupling value of the sum of couplings 
ri = -I- /ijv- The initial-values which lead to ri < (i.e., 
an ISL) become very small as N becomes large. We take both 
times t^/t = 0.95. 

B. Determination of the crossover energy 

Performing a numerical integration of the full set of RG 
equations (RGEs) , we calculate the crossover energy 
Ec — Ec{N), which separates ladder-like from 2D-AFM- 
like behavior (the RGEs are a generalization of the 3-leg 
RGEs given in Ref. the RGEs relevant for the AFM 
phase are given below). 

In Fig. H, we have plotted the strong-coupling value of 
the ratio of a 4-band AFM and a 2-band Cooper coupling, 
Ci2N-in/cin^ and of ri = c?i + /f^ for iV = 4, 8, 12, 16 
(we fix the strong-coupling value of at the band- 
width t). We take t±/t — 0.95, because for t±/t = 1, the 
band pair {l,N) is close to the van Hove singularities. 



such that its behavior would no more be typical for the 
system. When the number of chains, N, increases, the 
initial- value U/t which leads to a suppression of the 4- 
band-AFM couplings decreases fast. Similarly, for large 
N, the ri becomes only negative when U /t is very small. 
Note that for the other band pairs (j, j), the correspond- 
ing energy-scales are lower, since the velocities Vj are 
larger and the gaps are of the order of texp{—ajVj /U). 
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FIG. 3. The sum of couplings rj = c°j + fjj allows to de- 
termine between ISL and AFM: For < the band pair 
is an ISL and for > an AFM. We have plotted the 
initial-value Uc/t, which leads to ri ~ at the energy-scale, 
where the first coupling of the band pair (1, A'^) has grown up 
to the bandwidth t. The figure shows, that —logUc{N)/t is 
in good approximation a linear function of N (the squares are 
the calculated initial-values and the lines are a guide to the 
eye). 

Hence, both the 4-band-AFM processes and the rj 
have a strong and similar dependence on the energy-scale. 
For small enough initial- values U/t, the differences be- 
tween the Fermi velocities lead to a decoupling into band 
pairs. However, since Vj — vj+i cx 1/N, for increasing 
N, this decoupling becomes suppressed and the physics 
is dominated by 4-band-AFM processes. 

We calculate the crossover energy Ec as follows. We 
determine initial-values Uc = Uc{N) such that ri w 
at the scale where the first coupling of the band pair 
(l,iV) has grown up to the bandwidth t. For initial- 
values U < UciN), a spin-gap then opens, starting in 
the band pair (l,iV), and the 4-band-AFM processes 
become suppressed. We find that Uc{N) can be fit- 
ted by C/c(0)exp(-7iV), where 7 - 1 (for t^/t = 
0.8,0.85,0.9,0.95, we obtain 7 = 3.7,3.2,2.6,1.9) and 
Uc{0) ^ t, see Fig. |[ On the other hand, an initial-value 
U corresponds to an energy-scale (gap) t exp(~t/U). The 
crossover energy Ec is therefore given by 

Ec t exp[-a exp(7iV)] , (21) 
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where a ^ t/Uc{0) <C 1. Note that the energy Ec is 
also the upper Hmit for the spin-gap in the ISL phase. 
For large U (Heisenberg AFM), the corresponding scale 
(spin-gap) is J exp(— 0.68 A^), where J = AP/U (see also 
table I) g. While for large U the spin-gap decreases 
exponentially, in the small-C/ case the decrease is double- 
exponentially. 

As a result, the "phase" of the 2D system is for finite 
A^ present above an energy-scale Ec — Ec{N), where 
Ec{N) for N ^ oo. Equivalently, below the length- 
scale oc 1/ Ec, the dominating correlation functions are 
the same as in the 2D system. 



IV. THE 2D-LIKE AFM PHASE 

First (Sec. IV A), using the RG method, we calculate 
the asymptotic ratios of the couplings and the charge- 
gap of the phase present above Ec- Comparing our re- 
sults with the flow well away from half- filling jl^, we 
conclude, that AFM processes (at high energies) are nec- 
essary for an instability in the RG flow of the A^-leg 
Hubbard model. In Sec. IV B, we show that the phase 
above Ec is the quasi-lD counterpart of a 2D insulat- 
ing AFM. In particular, there is no odd-even effect and 
the interacting part of the Hamiltonian is similar to the 
2D Heisenberg AFM. Using bosonization techniques, we 
obtain, that the system is Mott insulating and that the 
insulator is of the same type as in the half-filled two- 
leg ladder, see Sec. IV C. For the spin-sector, we find 
the same physical properties as for the Heisenberg AFM, 
i.e., two gapless magnon modes, spinon confinement, and 
long-range order in the 2D limit. 



A. Asymptotic ratios of the couplings 

We give the RGEs of the 4-band-AFM couplings for 
large A^, i.e., keeping in the RGEs only sums over A^ 
products of couplings and dropping the (2-band) contri- 
butions of the order of 1. Wc then determine (analyti- 
cally) the asymptotic ratios of the couplings for Hubbard 
initial-values and the size of the charge-gap. 



1. RGEs 



N/2 

V - f c"- - 

y. y^jtlj kilk 



jkkj X ^ ^ t n n ^ 



dc' 



di ^ V, yj'^i'^ktjk iQ-m ktik 

N/2 



k ^ j 



jkkj 



dl 



dw 



jkk] 

dl 



X/ ( 2''^"-'"'="'= ' ^J^'^J^kizk ' '^kiik^jiij 
1 

~2 

X/ ~ ( ^^^mf^kilk + "jMj'^feMfc 



^yaf^kiik '^jrif^kiik ^kUk^'jiiJ 



(22) 



where the energy-scale is related to / by _E ^ te "^K 
Defining 



Jk "4^4=^ \jkkj '^jkkj 



we find from the above RGEs, that 

N/2 



Kk=j: H^4. 



(23) 



(24) 



i¥']-k 



Given the Hubbard initial- values 4c'', f_ = d^-.r- = 

jkkj Jkkj 

2uP^^-^_ = 2U/{N+1) and - 0, we obtain h+^{l) - 
for all I (this fixed point is stabilized by the 2-band in- 
teractions). Since 



dl 



i¥'j,k 



and h-f^{Q) > 0, the h-f^ flow to 0, h-j,{l) for increas- 
ing I. From the flow of the /i^^, we then calculate the 
flow of the 4-band couplings and obtain the asymptotic 
ratios 



t - 3g,k = 4c^,. = 3c 



L/.t- — iL. ., r_ — dc", r- = 4u'', r_— ^Su", r- (26) 
jkkj jkkj jkkj jkkj V ' 



TABLE I. The table shows the correspondence between 
crossover energy (spin-gap) and energy-scale of the 2D system 
for small and large U . The exponentially small charge-gap for 
small U leads to a double exponentially suppressed 1D-2D 
crossover energy (the result for the spin-gap in the large U 
case is from Ref. [6]). 



The RGEs of the 4-band interactions are given by a 
sum over A/2 products of 4-band couplings 



Jkkj _ y 1 ^ P 

L ^ Q,. V J^^3 kit, 



dl 



P P 
11 .. ir - 

k ' JiiJ kak 



i¥'j-k 

3 3 

"'"Te'^-'"-''^''"^ iQ^j^i^kiik 





2D energy-scale 


crossover energy 


u<.t 


texp{-\t/U) 


t exp[— a exp(77V)] 


u:>t 


J 


Jcxp(-0.68Af) 



6 



Note that u^^j.^, remains always small. The 2-band 
and single-band interactions are for large TV dominantly 
renormalized by the 4-band interactions. The asymptotic 
ratios of the 2-band and single-band couplings can there- 
fore be calculated by inserting the ratios ( p6| ) in the RGEs 
for these couplings . We then find, that the following 
couplings of the band pairs (j, j) grow and approach fixed 
ratios: 



^^3333 ^^3333 •^^3333- 



(27) 



The other 2-band and single-band couplings stay small; 
in particular the single-band SU(2) processes are small, 
|cjj | <C but they have become attractive, cj^ < 0. 

For bands k and j, which are close together on the 
FS, k j, the 4-band coupling becomes the same as the 
corresponding 2-band couplings, gjk ~ gj ~ gu- In the 
limit iV ^ oo, all the (2 and 4-band) ^-couplings take the 
same value, gj = .g,„ = g^i. The gap A is then similarly 
calculated as the asymptotic ratios; using that 



N/2 ^ 



i¥'3,k 



where 



Sjk 



3kk3 _ jfcfcj cf - u'' - 

4 4 'jfefej jfefej' 



we find for the scale of divergence 



In particular, — 2t/U for t±/t = and 

Tit 1 



In 



U ln(2iV/7r) 



(28) 



(29) 



(30) 



(31) 



for t ^ t±, where In < t/U for the validity of our calcu- 
lations. The logarithmic corrections come from the fact 
that vi — Vf^ ^ t/N (van Hove singularities). The gap 
becomes 



(32) 



where A is a function oi t±/t and of the order of 1 , A ~ 1 
(see also table I). 



2. Discussion 

It is instructive to compare our result with the situa- 
tion well away from half-filling, where the umklapp inter- 
actions can be neglected. In Ref. it has been shown, 
that without umklapp interactions and for repulsive in- 
teractions, U > 0, the large N limit is a FL, i.e., the 



typical gap-size scales exponentially to 0, oc . The 
reason is, that the forward scattering (f) processes, which 
tend to drive the system towards a RG instability give 
contributions of the order of 1, while the Cooper pro- 
cesses (c), which drive the system towards a FL, have a 
weight cx N. 

At half- filling, the situation is completely different: 
Here, the 4-band AFM processes with a weight cx N 
are (entirely) responsible for the RG instability at the 
energy-scale A ~ te~'^*/^ . Therefore, the A^-leg Hubbard 
model exhibits for large N only a RG instability, when 4- 
band AFM processes are (at least initially) present. The 
van Hove singularities, see Eq. (|3^), lead to a further 
increase of the gap-size (but are not the reason for the 
instability). 



B. Effective low-energy Hamiltonian 

We show that above Ec, the Hamiltonian of the half- 
filled A^-leg Hubbard ladders is similar to the 2D Heisen- 
berg AFM. 

At low energies (but above Ec), the ratios of the cou- 
plings are given by Eqs. ( p6| ) and (p7|), such that the U(l) 
and SU(2) 4-band interactions, see Eq. (^2|), simplify to 

~2^'^^^^,,,^l-^^^nks - (33) 
(for s =1, s =[ and vice versa) and 

- 4vl't^^.^vl;^^-,vl/t^^^_vl,^^^. (34) 

Defining 



n tak 

is 



-3 ^ Rj 

the Hamiltonian then takes the form 



,9i3 



1,3 



dx M, • Mo- 



(35) 



(36) 



The 2-band couplings /Jj'^ and u^^^^ give the contribu- 
tions Mj • Mj and the 2-band couplings c^i°^ and u''— 
lead to the products • Mj. The 4-band couplings 
result in the other products • Mj. For N large, all 
couplings take the same value gij = g > (s-wave AFM). 

As a result, all band pairs (j, j) are interacting with 
each other. In particular, for N odd, the band r = 
{N + l)/2 is interacting with all other band pairs and 
there is no qualitative difference between odd and even 
N. This contrasts the ladder-case at energies below Ec, 
where only interactions within the band pairs are 
present, respectively for odd N , within the band r. This 
then leads to an odd-even effect, i.e., the band r present 
only for odd N exhibits a gapless spinon-mode p9| . 

It is instructive to Fourier transform Eq. (pq) 
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H.c. 



(37) 



The Fourier transformed 2D AFM Heisenberg Hamilto- 
nian (large-?/ limit) takes the form 

Hj = -J J2 +6"^") 

k,k',q 

xM't^(k)rf^^,vI;,,(k+(7r,7r)+q) 
xvI/t^(kX,.'*4(k'-K'^)-q)' (38) 



where we substituted q by q-l-(7r, tt). Therefore, the inter- 
acting part of the Hubbard Hamiltonian ( |37| ) is basically 
the Heisenberg Hamiltonian restricted to the umklapp 
surface (the operators 'i'lfj^ and '^ljs have a momentum 
difference of tt) . Since U /t is small (and therefore the en- 
ergies low), there is no q-dependent coupling in Eq. (37). 
In real space, this corresponds to long-range spin-spin in- 
teractions. Note that the generalization/crossover of the 
Hamiltonian ( ^7|) to 2D is straightforward. Introducing 
a sum over qy is sufficient {q = qx)- 

To conclude this section: The half-filled A^-leg Hub- 
bard Hamiltonian — with on-site repulsion between elec- 
trons — becomes at energies below the gap A a Hamilto- 
nian with purely long-range AFM spin-spin interactions. 



C. Physical properties: Bosonization 

Using bosonization techniques, we derive the physical 
properties of the AFM Hamiltonian, Eq. (^). 



1. Bosonized Hamiltonian 



The noninteracting part is given by Hq, see Eq. (15). 
At the AFM fixed-point, the interacting part is as follows. 
The 4-band interactions take the form. 



H'^i = -9,k -He [msi +ns2) + (j ^j,k^k). (39) 



The 7^51 contains the spin-fields $0-77- and 0, 



Hsi = cos[/3($^jj_ - 



■fefe- 



kk- 



-)], (40) 



where f3 — and 7^52 contains the spin-fields ^an+ 



and Q. 



ns2 = cos[/3(6i<,jj_ - e^kk- 



- $^73+)]. (41) 



The charge-part. Tic, includes the charge-fields ^pjj+ 
and 9n 



He = cos[/3($pfcfc+ - 



pjj+ + "pjj 



"pkk-l\ 

'pn-~Gpk-k-)]- (42) 



In Tic, the first cosine comes from the 4-band non- 
umklapp interactions (p3|), while the second cosine re- 



sults from the 4-band umklapp interactions (34). The 



bosonized form of the 2-band non-umklapp terms is 



n 



2B 



-2cos(2/3$^jj+) - 4cos(2/3$„jj_) 
+2gj cos(2/3$o-jj-i-) cos(2/36'crjj-) 



(43) 



and of the 2-band umklapp terms 



^',2 = 9] cos(2/3$pjj+) [2cos(2/30^j,-_) 

-2cos(2^$<,jj+) - 4cos(2/?$^jj_)] 
-3gj cos(2/3$pjj+) cos(2/36'pjj_). 

The LLPs at the AFM fixed-point take the form 



(44) 



I TTV 



'3 T 3g,-/4 



TTV. 



± 35,/4 



(45) 



(implying Kpj+ < 1 and Kpj- > 1) and 



K, 



l 7rVj±gj/4 



(46) 



i.e., K^j+ > 1 and K^j- < 1. Table II gives an overview 
about the LLPs in the AFM and ISL phase. 

Note that the commutation relation for the field and 
its dual field hinders the pinning, i.e., the localization of 
the field and its dual field 9a in the minimum of a 
cosine at the same time. 

TABLE II. The phase at a particular energy-scale is re- 
lated to the value of the LLPs at this scale. The table shows 
the LLPs in the insulating AFM, ISL, and superconducting 
(SC) phase. Note that always ifpj- > 1 and K^j- < 1. The 
differences are, that Kpj+ < 1 for the AFM and ISL and 
Kaj+ < 1 for the ISL and the SC. 













AFM 


< 1 


> 1 


> 1 


< 1 


ISL 


< 1 


> 1 


< 1 


< 1 


SC 


> 1 


> 1 


< 1 


< 1 



•'PJJ-' 
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2. Charge-sector 

For the charge-sector, we then find, that the same 
fields are pinned as in the half-fiUed two-leg (respec- 
tively iV-leg) Hubbard ladder i.e., $pjj+ « and 
dpjj- ~ 0. The type of Mott insulator is thus the same in 
ID and 2D. Note that pinning of ^pjj- instead of Qpjj-- 
leads to another type of insulator. The difference between 
these two types of insulators becomes physically relevant 
upon (hole) doping. While the second type most likely 
becomes a FL, the first type becomes (in case we have a 
spin-gap) a superconductor, since doping only depins the 
$p-fields but does not destroy the phase coherence, i.e., 
the 0p-fields remain pinned, see Ref. []l9| . 

3. Spin-sector 

In contrast, in the spin-sector, both the fields $crjj± 
and the dual fields Ocrjj± appear in a cosine, resulting in 
a competition between different "phases", i.e., pinning of 
the field vs. its dual field. However, since Kcrj+ > 1 and 
K^j_ < 1 (for a comparison, see Sec. HI A), it is more 
favorable to pin O^jj^ and ^ajj- than the corresponding 
dual fields. 

The physical interpretation is then the following. The 
pinning of ^ajj- leads to a spinon confinement, leaving 
as physical particles spin 1 magnons. Note that only the 
differences 9cjj+ — Qakk+ appear in H'^^ , such that the 
total magnon mode(s), given by 9t = \/'2'/N dcrjj+ 

and <I>T = ^crjj-i-^ remains gapless. The n and 

mode then result from a superposition of left/right go- 
ing modes along the chains (since we use open boundary 
conditions, the transvers momentum is always positive). 

For the calculation of the spin-spin correlation func- 
tion, we first rewrite the real-space operators dig in terms 
of the band operators '^hjs and then the band operators 
in terms of the bose-operators. Using that only products 
of fields, which contain the pinned charge-fields ^pjj+ 
and 9pjj- give non-vanishing contributions to the spin- 
spin correlation function, the real-space spin-operator at 
the position [x, i) takes the form 

j ^ 

+i? ^ L ) , (47) 



where jjm = ■\/2/(A^ + 1) sm[TTjm/{N + 1)] and at 
half- filling, kpj + kpj — tt. The abelian bosonization 
scheme used here, breaks non-abelian symmetries, i.e., 
the SU(2) spin-symmetry is broken down to U(l), see, 
e.g., Refs. 1^,^. In our case, only the x and y com- 
ponents of the spin-spin correlation function give then 
straightforwardly the correct physical result (this "prob- 
lem" occurs also in a single chain with a spin-gap [p2t). 



The products of fields appearing in Eq. ( ^7| ) are then 
rewritten in terms of bose-operators according to 

'^Rj^'^Lji oc e'^("*''"+~*-""-+^«^"-+^''"+^ (48) 

The charge-fields and ^ajj- are pinned and can be set to 
zero. Using that 0a-j]+ — ^'cr/cfc+ ~ 0, we express Oajj-t- in 
terms of the total spin- mode 9t, Ocrjj+ = y/2/N9T- Since 
Hq is gaussian in the fields and 9t, the spin-spin cor- 
relation function takes the forrn (for details about such 
calculations, consult, e.g., Ref. pTj ) 

{S^{x) ■ Si{0)) OC (-1)'+' cos(7rx)/xi/^, (49) 

where x < oc 1/Ec (for x > we are in the ISL phase, 
where the decay is exponentially with coherence length 
S,AF S,ci respectively, for N odd, cx 1/x). 

The energy-scales for charge and spin-excitations are 
different. The LLPs in the charge-sector, see Eq. (^), 
deviate more from their noninteracting value {— 1) than 
the LLPs in the spin-sector, Eq. (^), i.e., 

\Kpj± - 1| ~ igj/i^^Vj), \K,j± - 1| « gj/i^irv,). (50) 

This implies that the charge-fields are pinned more 
strongly than the spin-fields and therefore the Mott-gap 
(we may interpret the RG-scale A as the charge-gap) is 
larger than the gap to spin-1/2 (spinon) excitations. The 
coupling between the different bands makes it difficult to 
investigate the detailed excitation spectrum, such that 
we do not go beyond this qualitative result (for the pure 
sine-Gordon model, there exist rigorous results about the 
gap dependence on Ka, see Ref. p7[|). 

V. EFFECT OF (HOLE) DOPING 

Next, we study the effect of doping away from half- 
filling. There are two different doping regimes. The 
lightly doped case, where the influence of the umklapp 
processes has to be taken into account, can be treated as 
a perturbation of the half-filled low-energy Hamiltonian 
(limit (5^0 and U finite). For increasing doping, one 
has then again to investigate the RG flow. 

For the lightly doped case (Sec. V A), we find above Ec 
the same phase as below (see Ref. |Q ) : a conducting 
spin-liquid (bound hole-pairs), see also Fig. |^. Further- 
more, while at half-filling the magnetic correlations are 
only at (tt, tt), they become incommensurate upon dop- 
ing. 

In Sec. V B, we briefly show how and why for increas- 
ing doping phase-coherence between all band pairs sets 
in and the system becomes a 2D-like d-wave supercon- 
ductor. The low-energy Hamiltonian takes the form of 
a d-wave BCS model, where the pairing is between elec- 
trons. 

Note that the superconducting instability is for N ^ 2 
generated in the same way as for the 2D case [|4|-^: 
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A Kohn-Luttinger-type attraction [g^ is mediated by 
AFM processes In contrast, the superconducting 

instability in the two -leg ladder is generated by forward- 
scattering processes Q. 



A. Lightly doped case 



The ISL phase of the half-filled two-leg Hubbard lad- 
der becomes superconducting when doped, because of the 
phase-coherence (see discussion of the charge-sector in 
the previous section) and the spin-gap present at half- 
filling |^6|. Here, we show that the doped Hubbard 
AFM becomes a conducting spin-liquid, formed by bound 
hole-pairs (this phase may becomes superconducting at 
very low energies; for the 4-leg ladder, see discussion in 
Ref. @). 

The 2-hand interactions alone, see Eqs. (^ ) and (44), 
would result in an ISL phase, i.e., '^ajj± ~ 0. It is the 



presence of interactions between bands i and j which are 
not "paired" (j ^ which renders the system a 2D- 
like AFM. The effect of a decreasing TV < oo is therefore 
that interactions between unpaired bands become sup- 
pressed and the AFM correlations are reduced. Next, we 
show that the effect of doping away from half-filling is a 
similar one, i.e., the doping i5 corresponds to 

At low dopings and at finite N, we can dope the system 
perturbatively by introducing a chemical potential term 
— /iQ (Q is the total charge) in the half-filled low-energy 
Hamiltonian, see Ref. For the charge-sector, we can 
use bosonization, where the charge density is given by 




71 71 

FIG. 4. The square is the umklapp surface, which is 
the FS at half- filling (t — t±). Left; There are 2 types of 
AFM processes; AFM processes which are umklapp processes 
(dashed arrows), and AFM processes which are not umklapp 
processes (solid arrows). Right: The (non-umklapp) AFM 
processes which take place within a band pair {j, j) are iden- 
tical to the corresponding (d-wave) Cooper processes within 
this band pair. Upon doping away from half-filling, these 
processes remain large, open a spin-gap and lead to phase 
coherence within band pairs (j,j). Note that we use for the 
ladders open boundary conditions, whereas the above figures 
are drawn — for a better comparison with the 2D case — for 
periodic boundary conditions. 



Doping then introduces kinks in the fields ^pj]+ , such 
that the expectation value of cosine terms, which contain 
this field gradually goes to zero (the kink is a fermionic, 
strongly localized particle). Here, the kinks are bound 
hole-pairs of charge 2 and zero spin. All 4-band interac- 
tions contain the field ^pjj+, see Eq. (p^), and therefore 
vanish. The effect of the doping on the 2-band inter- 
actions is the same as when doping an ISL, i.e., in the 
bosonized version, the umklapp term (Q) vanishes leav- 
ing only the term (^3|), which opens a spin-gap, and leads 
to d- wave-like phase coherence between the bands j and 
J, corresponding to the Fermi momenta {kpj, n—kFj) and 
(tt — kpjjkpj)- The charge-gaps close to (7r/2,7r/2) are 
the smallest ones (note that gN/2/vN/2 < ■■■ < gi/vi), 
such that the hole-pairs enter first there. 

In other words, it is the fact, that the (s-wave) AFM 
processes 



*fl,i*iJT*LjT*«J"i (+other terms) 



(51) 



(partially) coincide with the (c?-wave) Cooper processes 
with momentum transfer tt 

^l^J = *ki*LT*iJT*flJi (+othcr terms), (52) 

which makes a doped AFM a superconductor: These par- 
ticular AFM/ Cooper processes flow at higher energies 
— driven by AFM processes — to strong coupling and 
are not suppressed at lower energies by doping (see also 
Fig. U). We therefore have AFM mediated superconduc- 
tivity |24| . The full phase-coherence around the FS then 
grows out of the phase-coherence between bands j and j, 
see below. 

Upon doping, the AFM correlations become mcom- 
mensurate. When doping away from half-filling, kpj + 
kpj < TT, and the AFM peak shifts away from (tt, tt), see 
Eq. (^7|). In particular, for small tj^/t <C 1, all bands are 
(almost) equally doped and kpj + kpj = 7r(l — 5), i.e., 
the incommensurability is equal the doping 5. Further- 
more, doping leads to a cutoff of the AFM processes at 
the energy-scale of the chemical potential /i and therefore 
to an exponential decay of the spin-spin correlation func- 
tion, where the coherence-length is given by ^af oc l//i. 



B. Increasing doping 

We argue, that for increasing doping, phase-coherence 
between the band pairs sets in, when the "distance" 
between neighboring bands is of the order of the "dis- 
tance" of the FS to the umklapp surface, i.e., for dopings 
5 > Sc{N) ~ {t±/t)/N. The system then becomes a 
2D-like c?-wave superconductor. 

Since a finite chemical potential fj, results in a low- 
energy cutoff for the umklapp and 4-band- AFM interac- 
tions, one has to start initially with all interactions which 
are large at this cutoff energy-scale. In leading order, this 
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implies that not only the 4-band interactions (f-fj.- and 

JKKJ 

^jkkj '^^^^ to be taken into account, but also "neighbor" 
processes of the form c^^'^_|_^__|_^ and These 
processes are only present down to an energy-scale tj_/N. 
For dopings < 6 < {t±/t)/N, the "neighbor" processes 
are cutoff before the main AFM processes, which leads to 
a decoupling into band pairs {j, j) , while for 5 > (t /t) /N 
all AFM processes are cutoff by the chemical potential. 
In this case, not only 2-band Cooper interactions within 
band pairs (j, j) are large at the cutoff energy-scale, but 
also the neighboring 2-band Cooper interactions (j, j±l). 

Integrating the RGEs without AFM and uniklapp pro- 
cesses, but with initial- values for the Cooper and forward 
interactions given by the AFM phase, see Eq. (p7|), and 
similar initial-values for non-leading Cooper processes 
Cj'j'±p (j> = 1,2,3,.. .), we obtain a dx2-y2 superconduct- 
ing instability at lower energies. We always find that the 
scaling of the couplings is exactly such that (in particu- 
lar, 44=c^,/(;«/5«0) 

H^Ho+J2 [ '^^ ' (53) 

where the pairing operator is Aj = ^'_Rj| ^l^x ^ Rji 

and Vij < for i,j < N/2 respectively i,j > N/2 and 
Vij > in all other cases {d^2_y2 symmetry). The size 
(but not the sign) of Vij depends on the initial-values 
(i.e., on the doping). Note that the energy-scale of the 
system is reduced by doping, i.e., the superconductor has 
a lower scale than the AFM. 

For large dopings, the 4-band interactions are cutoff be- 
fore they could have flown sufficiently close to the AFM 
fixed-point, such that there is no instability in the RG 
fiow and the low-energy phase is a FL ||l^ . 

As a result, for N finite, there is between the insulating 
AFM respectively ISL phase at half-filling and the d-wave 
superconductor at higher doping, an intermediate, con- 
ducting phase with a spin-gap (see also Fig. |l|). In the 
2D-limit, the intermediate phase seems to disappear. 



VI. CONCLUSIONS 

We have studied the dimensional crossover from ID 
to 2D in half-filled, weakly interacting, A^-leg Hubbard 
ladders for t± < t. Using the RG method, we found a 
crossover energy Ec texpl—aexp^-^N)], where at ener- 
gies below Ec the behavior of the spin-sector is ladder-like 
and at energies above Ec 2D-AFM-like. We showed that 
the 4-band- AFM processes are responsible for an instabil- 
ity in the RG flow. Since without 4-band and umklapp 
processes the Hubbard ladders become a FL for large 
N [Q, AFM processes (at higher energies) are necessary 
and sufficient for a nontrivial phase in the (doped) TV-leg 
Hubbard model. 

In the AFM phase, the couplings flow towards uni- 
versal ratios and we obtained an analytical expression 



for the Hamiltonian which is similar to the Heisenberg 
AFM. Using bosonization techniques, we further investi- 
gated this Hamiltonian and found that the charge-sector 
is the same as in the half-filled two- leg ladder, i.e., the 
type of Mott insulator is the same in ID and 2D. For 
the spin-sector, bosonization techniques confirmed that 
the small-[/ case has similar physical properties as the 
large-[/ (Heisenberg) case. 

The effects of doping on the half-filled Hubbard Hamil- 
tonian and of reducing the number of coupled chains N 
is (almost) the same: In both cases, the interactions 
between unpaired bands are suppressed, the antiferro- 
magnetic correlations are reduced and only interactions 
within band pairs remain. Doping then opens a 

spin-gap, but phase coherence exists only between bands 
j and J. At higher doping levels 5 ~ {t±/t)/N, we can 
expect that phase coherence between all bands sets in, 
such that the system becomes a 2D-like d-wave super- 
conductor. While the lightly doped case is conveniently 
described by bound hole-pairs, the 2D-like superconduc- 
tor at higher dopings is represented by a BCS-like d-wave 
Hamiltonian with electron pairing. In the doped case the 
antiferromagnetic correlations become incommensurate. 

As a result, the Hubbard model cannot be under- 
stood in terms of an AFM, ISL or superconductor alone. 
The phase (correlation function) depends always on the 
energy-scale (respectively length-scale) and is related to 
the value of the LLPs. The AFM exhibits the highest 
energy-scale. At half-filling, the groundstate (for N fi- 
nite) is the ISL and for sufficient doping the d-wave su- 
perconductor. We conjecture that in the doped case, the 
fiow of the couplings from the AFM to the d-wave takes 
place via an intermediate phase: Consult table II; AFM 
and superconductor differ in two LLPs. This intermedi- 
ate (finite energy) phase has then the same LLPs as the 
ISL — suggesting that the spin-gap opens before phase 
coherence sets in. 

We believe, that also for other systems of interacting 
electrons the dimensional crossover from ID to 2D can 
be studied in a similar way and maybe enlightens our 
understanding of the 2D case. 
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